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Abstract — Image fusion refers the process of combining 
multiple images of a scene to obtain a single composite image. 
The composite image should contain a more useful description of 
the scene than provided by any of the individual source images. 
Image fusion is widely used in intelligent robots, stereo camera 
fusion, medical imaging, and manufacture process monitoring, 
electronic circuit design and inspection, complex machine/device 
diagnostics and in intelligent robots on assembly lines. 

Here we have consider some image fusion techniques such as 
averaging, min-max, PCA, brovey, pyramid based and 
transform based techniques. Various parameters have been 
considered and varied for all these techniques and found that 
every algorithm has some advantages and drawbacks. 
Combination of qualitative and quantitative assessment 
approach may be the correct way to find out which fusion 
algorithm is most appropriate for an application. 

Index Terms — Image Fusion, Spatial Domain, Frequency 
Domain Techniques, Wavelets etc. 


I. Introduction 

Latest technologies in image capturing devices help us to 
extract variety of different information from an image. This 
information can be collectively combined using “fusion” to 
generate more informative image. Multi-view, multi-modal, 
multi-temporal and multi-focus are the four ways in which 
image fusion can be performed [1]. Image Fusion reads 
several images of the same scene or objects and retrieves 
important information from them to put it into a single output 
image. The basic idea is to perform a multiscale transform 
(MST) on each source image, then construct a composite 
multiscale representation from these according to some 
specific fusion rules [2]. 



There are various image fusion techniques which are 
classified below. 



Image Fusion Techniques 

Averaging Technique: 

It is a well documented fact that regions of images that are in 
focus tend to be of higher pixel intensity. Thus this algorithm 
is a simple way of obtaining an output image with all regions 
in focus. The value of the pixel P (i, j) of each image is taken 
and added [3]. This sum is then divided by 2 to obtain the 
average. The average value is assigned to the corresponding 
pixel of the output image which is given in equation (2). This 
is repeated for all pixel values. 

The fused image K( i,j) is given as 

K (i, j) = {X (i, j) + Y (i, j)}/2 (1) 

Where X (i, j) and Y (i, j) are two input images and K(i,j) is 
the fused image. 

Maximum Selection Scheme: 


Image fusion scheme This sc heme just picks coefficient in each subband with 

largest magnitude. A selection process is performed here 
wherein, for every corresponding pixel in the input images, 
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the pixel with maximum intensity is selected, and is put in as 
the resultant pixel of the fused image 

K( i,j) K(i, j) = Max.[w(H(x,y)), w(I2(x,y))] (2) 

Where Il(x,y), I2(x,y) are the input images. 

Minimum Selection Scheme: 

This scheme just picks coefficient in each subband with 
smallest magnitude. A selection process is performed here 
wherein, for every corresponding pixel in the input images, 
the pixel with minimum intensity is selected and is put in as 
the resultant pixel of the fused image [4]. 

K( i,j) K( i, j) = Min.[w(Il(x,y)), w(I2(x,y))] (3) 

% Load two original images: a mask and a bust 
load mask; XI = X; 
load bust; X2 = X; 

% Merge the two images from wavelet decompositions at 
level 5 

% using db2 by taking two different fusion methods 

% fusion by taking the mean for both approximations and 
details 

XFUSmean = wfusimg(Xl,X2,'db2',5,'meanVmean'); 

% fusion by taking the maximum for approximations and the 
% minimum for the details 

XFUSmaxmin = wfusimg(Xl,X2,'db2',5,'max','min'); 

Principal Component Analysis: 

It is a simple non-parametric method of extracting relevant 
information from confusing data sets. PC A is a useful 
statistical technique that has found application in fields such 
as face recognition and image compression, and is a common 
technique for finding patterns in data of high dimension [5]. 
The origins of PC A lie in multivariate data analysis, it has a 
wide range of other applications PC A has been called, ’one of 
the most important results from applied linear algebra and 
perhaps its most common use is as the first step in trying to 
analyses large data sets. 



PCA process flow 


The input images Ii(x, y) and I 2 (x, y) are arranged in two 
column vectors and their empirical means are subtracted. The 
resulting vector has a dimension of n x 2, where n is length of 
the each image vector. Compute the eigenvector and 
eigenvalues for this resulting vector are computed and the 
eigenvectors corresponding to the larger eigenvalue obtained 
[6]. The normalized components Pi and P 2 are computed from 
the obtained eigenvector. The fused image is given by 
equation, 

I f (x,y)=P 1 Ii(x,y)+P 2 I 2 (x,y) (4) 

coeff = pca(X) 
coeff = pca(X,Name, Value) 

[coeff, score, latent] = pca(X) 

[coeff, score, latent, tsquared] = pca(X) 

[coeff,score,latent,tsquared,explained,mu] = pca(X) 

Brovey: 

It is also known as the color normalization transform because 
it involves a red-green-blue (RGB) color transform method. 
The Brovey transformation was developed to avoid the 
disadvantages of the multiplicative method. It is a simple 
method for combining data from different sensors. It is a 
combination of arithmetic operations and normalizes the 
spectral bands before they are multiplied with the 
panchromatic image [7]. It retains the corresponding spectral 
feature of each pixel, and transforms all the luminance 
information into a panchromatic image of high resolution. 

Red = (bandl/S band n)* High Resolution Band 

Green = (band2/£ band n)* High Resolution Band 

Blue = (band3/£ band n)* High Resolution Band 

High resolution band = PAN (5) 

Pyramid: 

A pyramid based fusion comprises of a number of images at 
different scales which together represent the original image. 
Every pyramid based image fusion has three stages. Number 
of levels ‘L’ of pyramid is pre decided based on the size of the 
image. 

The image reduction process involves lowpass filtering and 
downsampling the image pixels. The image expansion 
process involves upsampling the image pixels and lowpass 
filtering. You can also use this block to build a Laplacian 
pyramid. 

A pyramid structure contains different levels of an original 
image. These levels are obtained recursively by filtering the 
lower level image with a low-pass filter [8]. We first make a 
Gaussian pyramid by filtering each level of image using a low 
pass filter and do the down sampling. As the level goes up, the 
image is getting smaller and smaller. 

The matlab code is as follows: 

function LLk = Lk(k,address) 
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Now, we do the reconstruction part to reconstruct the original 
image using the first level of Laplacian pyramid and the 
filtered, upsampled version of (k+l)-th level of Gaussian 
pyramid [9]. 
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Reducel Expaixll Laplacian (bi£1) 

Pyramid Flow chart 
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Conclusion: 


I have compared various image fusion method based on a 
Gaussian mixture distortion model. The results showed the 
advantages of pyramid approach in some cases. We also have 
studied the effect of the settings (number of Laplacian 
pyramid levels, the size of local analysis window). Generally 
speaking, using more Laplacian pyramid levels can be 
beneficial but this comes at the cost of increased complexity. 
In practice, we have found there is not much difference 
between the fused results obtained using a different number of 
pyramid levels if a sufficient number of levels is used (5 in our 
experiments). The local analysis window should be small 
enough so that the parameters are indeed constant in the 
window, but it should be large enough to contain enough 
sensor data to estimate the parameters reliably. In cases we 
considered, we found a 5x5 or 3x3 window size is a good 
choice. 
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